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Abstract 

We calculate the bubble nucleation rate in a fourth-order gravity theory 
whose action contains an R 2 term, under the thin-wall approximation, by 
two different methods. First the bounce solution is found for the Euclidean 
version of the original action. Next, we use a conformal transformation 
to transform the theory into general relativity with an additional scalar 
field and then proceed to find the bounce solution. The same results are 
obtained from both calculations. 



1. Introduction 

Consider the theory of a single scalar field defined by the action 



S = / <fx 



i<V>^0-\/(0) 



'1.1} 



where V has two local minima, cf>±. One of these, 0_, is a global minimum ("true 
vacuum"). The other corresponds to a metastable state ("false vacuum") that 
decays through bubble nucleation. The decay rate from the false vacuum to the 
true vacuum per unit volume can be calculated in the semiclassical limit 

r/v = A e - B ^[i + o(h)]. (1.2) 

Algorithms for computing the coefficients B and A have been given by Coleman[l] 
and Callan and Coleman [2] respectively. 

In applications to the early universe, the effect of gravity becomes interesting. 
When gravity is weak, the basic features of bubble nucleation are similar to those 
in flat space, with small modification of the decay rate and the bubble radius. If 
we confine ourselves to the case where both the true and false vacua have positive 
energy density, then the true and false vacua will be de Sitter spaces. Gravita- 
tion makes the materialization of the bubble more likely, and makes the radius of 
the bubble at its moment of materialization smaller. Coleman and De Luccia[3] 
have calculated the coefficient B in the case where gravity is described by general 
relativity. 

The Hilbert action for general relativity contains a single term proportional to 
the curvature scalar, and leads to second-order field equations. For some modified 
theories of gravity, the action contains terms that are combinations of the curvature 
scalar, Ricci tensor, Weyl tensor and the full Riemann curvature tensor. This 
leads to higher-order field equations. Such theories have been considered (see, e.g., 
[7]), for example, in the context of variations of extended inflation [8, 9]. In this 



paper, we will confine ourselves to the case where the gravitational action contains 
an arbitrary function of the curvature scalar. The method can be easily applied 
to the case where the action contains arbitrary functions of the curvature scalar 
multiplied by a function of another scalar field. 

When there is only an additional R 2 term in the action, Whitt found a confor- 
mal transformation that transforms the theory to standard general relativity but 
introduces a new scalar field [5]. This method can be generalized to the case of an 
arbitrary function of the curvature scalar. This is done in Section II, where some 
general features of fourth-order gravity are discussed. To calculate the bubble nu- 
cleation rate in fourth-order gravity theory, one could use the Euclidean version of 
the original action. This action contains higher powers of the curvature scalar. The 
bounce solution[l] of this Euclidean action gives the coefficient B of the bubble 
nucleation rate. Alternatively one could use a conformal transformation to trans- 
form the theory into general relativity with an additional scalar field. One then 
could proceed following Coleman and De Luccia[3]. In Section III, we calculate 
the coefficient B when the extra terms are small, using a perturbation method. In 
Section IV, we will use the conformal transformation to calculate the coefficient B. 
The result of this approach is found to agree exactly with the direct approach. We 
have made no attempt to calculate the coefficient A. Even in general relativity, 
there is not a satisfactory method to evaluate this coefficient. 



2. Fourth-Order Gravity and Conformal Transformation 

Consider a theory described by the action 



,4 ,—_ I W) 



where R is the curvature scalar and C m is the Lagrangian of the matter fields. 
General relativity is the case when f(R) = R. For simplicity, we will confine 
ourselves to only one matter field with the matter Lagrangian 



C m = -d lt (f>d»<!>-V((i>) (2.2). 

The generalization to more fields is straightforward. It is easy to show that the 
equations of motion for this theory are 

\f{R)g^ - f(R)Rp, + f"(R)(R^ - iVV) (2 3) 

+ f'"(R)(R ;fl R. u - R ;a R ;a 9^) + 8ttG(^0 d v <t> - C m g^) = 



and 

^;A ;A = -^ir^- (2-4) 



; a dV{4>) 



Because the curvature contains second order derivatives of the metric, the above 
equations contain fourth-order derivatives of the metric unless f"(R) = or, in 
other words, if the theory is general relativity. However, the higher order terms 
enter the equation only through the derivatives of the curvature scalar. We can 
lower the order of the differential equations by two if we introduce a new variable 
a and arrange for this to be equal to the curvature scalar. The equations are then 
second order, but we have to supply the equation 

a(x) = R(g^) (2.5) 

for the new variable. 



This can also be done by looking at the action directly. The action (2.1) is 
equivalent to 



S 



„/4 / 

a xy—g 



!6nG 



[f'(a)R-af'(a) + f(a)}+jC r 



By varying a we get 



<r = R 



(2.6) 



(2.7). 



Substituting this relation back into the action (2.6), we find that it is exactly the 
same as the original action (2.1). To make the theory look like general relativity, 
we introduce the conformal transformation 



9ixv = f'(cr)g^. 



(2.8) 



In the new metric the connection is 



■pA "pA i y^A 



(2.9) 



where the tensor Ey^ is 



.A 



/» 



0,n9v + v,v9^ - v ,X 9vv 



»" If [a 
The Ricci tensor in the new metric is found to be 



-H-flU — -tiflV 



V'{v) 



X 



2o>, + o- ;A fiV 



r 

v 



1<r,H0,v + 0",AO" 9n» 



3/7" 
2 V/' 



+ - I — I 0>C> 



and the new curvature scalar is 



r =R^r v 

1 



R 



sr ,a , afr\ 2 jA 



^ a ' Xx -— a > Xa '" + 2{j 



(2.10) 



(2.11) 



(2.12) 



where the indices are raised and lowered by the old metric. Also, 

^ ; \|old = o,^ v = (fa>\ - f"a, x a> X )\ new 
C,AC' I old — f a ,X a ' Inew 

With the aid of the above relation, the action can be written as 

R 1 



(2.13) 



S = I d A x 



16nG 16nG 



3 ff'V x cr f 

7) °"* + r-F 



d^dycjyg^ V((l>)\ 

V f 2 J 



(2-14) 



where the arguments of /, /' and /" are a. Now we have the action of general 
relativity coupled to two scalar fields, although the matter Lagrangian is somewhat 
unconventional. It can be shown that the equations of motion for the new action 
are equivalent to those of the original action, which proves the legitimacy of the 
above procedure. 

Before we go to the next section, let us look at the de Sitter space solution of 
fourth-order gravity. De Sitter space is a homogeneous space with 0(4, 1) symme- 
try. Because of this symmetry, the Ricci tensor can be written as 



iV = l R 9^ (2-15) 



where R is a constant throughout the whole space-time. The matter field is also 
homogeneous throughout the space-time manifold, so the equation of motion (2.3) 
reduces to 

\f{R)R-\f{R)=KV{<h) (2.16) 

where 0o is a stationary point of V and k = 8ttG. This is an algebraic equation 
which can be solved for R. 



Let us consider the specific case where the action is 

R 



S = J d*xV=9 [~Y^g + aR2 + 2 9 ^ 9 ^ ~ V{(P) \- (2,17) 

Then /(a) is 

f(a) = a - !6nGaa 2 , (2.18) 

and the conformal transformation is 

9ixv = (1 - AaKa)g^ v (2.19) 

In this new metric the action is 

. , R | 12a 2 Kfy<rfflV <V><^0 V(0) + mx 



2k (l-4a/to-) 2 2(1 -4ana) (1 - Aana) 2 

'(2.20) 

For both general relativity and the theory described by (2.17), equation (2.16) 
gives 

R = -4/d/(0 o ) (2.21) 

because the a terms happen to cancel. If higher order terms enter the action, the 

result differs in general from the general relativistic result. In the case when only 

R and R 2 terms enter, this solution corresponds to the stationary point of the last 

term of (2.20), 

V{(j>o) + aa 2 
(1 — 4aK,a) z 

This stationary point is located at a — — 4nV((j)o), which is just the curvature 
scalar in the original metric. The curvature scalar in the new metric is 

R= 4 " V y, , (2.23) 

in agreement with Eq. (2.12). 



3. Bounce Action: The Direct Approach 

Consider a theory described by the action (2.1), where V(<f>) has two local 
minima, (f>±, only one of which, 0_, is an absolute minimum. Let us also assume 
that both V ((/>-) and V((f>+) are positive. The classical field theory defined by this 
action possesses two stable homogeneous equilibrium states, = 0+ and = 0_. 
In the quantum version of the theory, though, only the second one corresponds to a 
truly stable state, a true vacuum. The first decays through barrier penetration; it 
is a false vacuum. The decay rate per unit time per unit volume can be calculated 
in the semiclassical approximation as 

T/V = Ae- B / h (l + 0(k)). (3.1) 

To calculate the coefficient B, let us consider the Euclidean version of the 
theory. The Euclidean action is defined as minus the formal analytic continuation 
to imaginary time of the Lorentzian action (2.17) 

s E = J d ' x V~g{- f -^ + \d^ M gT + v{4>) } (3.2) 

where the metric is the usual positive-definite one of Euclidean four-space. Let 
(4>,g) be a solution of the Euler-Lagrange equations associated with S E such that: 
[i] (0, g) approaches the false vacuum solution at large Euclidean distance, [ii] (0, g) 
is not identical to the false vacuum solution, and [iii] (0, g) has Euclidean action 
less than or equal to that of any other solution obeying [i] and [ii]. Then the 
coefficient B in the vacuum decay amplitude is given by 

B = S E (<f>,g)-S E (<i>+,g+) (3.3) 

where (<p + ,g + ) is the false vacuum solution, ((f), g) is called the bounce solution. 



It can be shown that in flat space the bounce is always 0(4)-symmetric[4]. 
In curved space, although not proven, it is very plausible that the bounce is still 
0(4)-symmetric. We will work under this assumption. The metric can then be 
written as 



ds 2 = d£ 2 + p 2 {£)d& 



(3.4) 



where £ is a radial coordinate from the center of the 0(4)-symmetry and p(£) 
measures the circumference divided by 2n at radial coordinate £. 

In the following, we will confine ourselves to the case where only R and R 2 
terms enter the action and the theory is described by the action (2.17). In 0(4)- 
symmetric coordinates (3.4), the equations of motion are 



,2 



1 + 



Kp (V- 

3 \2 V 


-V 






+ 36a 


"l + 2p' 2 - 

P 4 


-3p 1p p 

P 6 


p" 2 , y P " r 



(3.5) 



and 



„ 3p' , ^ 

p a0 



(3.6) 



We will use perturbative methods to find the effect of small a. To be self- 
consistent, the second term in the bracket should be much smaller than the first 
term. This is always true if 

9 
m 

a Ml <<l 

because the second term is at most of the order am® /M 2 where m is the mass 
scale of the potential V(4>). 

Under this condition, the Euclidean action can be split into an unperturbed 



part 

S = J d 4 x^g |~ + l -{d(frf + V(0)} (3.7) 

and a perturbation 

S a = I d A x^g{-aR 2 ). (3.8) 

Let (po,0o) be the bounce solution when a vanishes and (pi,0i) be a small 

correction 

P = po + pi 

(3.9) 

= 00 + 01 

Because (po, 0o) is a solution to the Euler-Lagrange equations of Sq, the first order 
correction to So due to (pi,0i) vanishes, so to first order the Euclidean action is 

Se = S [po,H + S a [po,H + 0(a 2 ) (3.10) 

and the bounce action is 

B = Bq + B a 

= {So(po,(f)o) — Sb(false vacuum)} + {S a (po, (fro) — ^(false vacuum)} 

(3.11) 

We will use the thin wall approximation in the following calculation [1,3]. The 
thin wall approximation is valid when the energy density difference between the 
false and true vacua is small compared to the potential barrier. In this case the 
thickness of the wall is small compared to the radius of the bubble. We will also 
assume that gravity is weak, i.e., that the scale m of V((fr) is much smaller than the 
Planck mass. The Euclidean true and false vacua with positive cosmological terms 
are Euclidean de Sitter spaces, which are four-spheres. We will further assume that 
the radii of these four-spheres are much larger than the flat-space bubble radius. 
Under above assumptions the bounce solution has a thin 3-dimensional spherical 
wall which separates an interior true vacuum region and an exterior false vacuum 



region. The radius of this wall is large compared to the wall thickness, but much 
smaller than the radii of the true and false vacuum 4-spheres. 

First let us consider the unperturbed solution with a — 0. Inside (outside) the 
bubble (po, 0o) is just de Sitter space of true (false) vacuum 

00 = 0± 

sin(if ± (3.12) 

with 

R± = Y1E\ = AkV{4>±). (3.13) 

This is the exact result for the true and false vacuum 4-spheres even when a ^ 
because the a terms happen to cancel. The contributions to Bq have been 
calculated by Coleman and De Luccia[3] 



£o| out = 

3 / 2 ] (3.14) 



UlT 2 



K 2 V((j)-) 



i - \kpIv{4>-) 



where 2-7rpo is the circumference of the 3-sphere thin wall. In the thin wall approx- 
imation, the contribution to Bq from the shell can be shown to be [3] 



5 | wall = An 2 pl J d^V^) - V(<j> + )] = 2tt 2 p^i 
where po is determined by requiring that the bounce action 

B = B o\ in + 5o| wa ll + A) lout 

be stationary. 



(3.15) 
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We need the curvature scalar in the wall region for latter calculations. Because 
(PO) 00 ) is a solution to the theory defined by So (which is standard Euclidean 
general relativity), (po,0o) satisfies Einstein's equation 



Rfiu ~ ^9^R = -kT^ v (3.16) 



The trace of this, 



gives 



^A 



R = kT a x , (3.17) 



R = K 



o 2 + 4V(0o) • (3.18) 



Because the difference between V(<p+) and V((f)-) is very small compared to V{4>) 
in the thin-wall region, define 

Vo = \[v(M + v^)] (3ig) 

[/(</.) = v(<j>) - v 

Except in the thin wall region, is almost a constant (either 0+ or </>_) and V is well 
approximated by Vq. Because the bubble's circumference radius is large, we can 
neglect the -%-<f>' term in Eq. (3.6). We then have the "approximate conservation 
law" 

i0 o 2 - V(<f> ) = -V . (3.20) 



We find 



o = v / 2W o 7 

R = k[6U(<I>o)+Wo\. 



(3.21) 



Now consider the first order correction to the coefficient B 

B a = S a (po, 0o) — S a (false vacuum) (3.22) 

In the exterior of the bubble the bounce solution coincides with the false vacuum 

11 



solution, so the contribution to B a vanishes 



5 "lout = °- 



(3.23) 



The contribution of the wall region to B a is 



B, 



a I wall 



—a 



d X [R - -Rfaise) 



(3.24) 



Using (3.21) we find that 



Taiwan = (-«k)2tt 2 ^ J ' d£ 6U(6U + 8Vo) 

= (-aK)2ir 2 pl f -jL QU(QU + 8Vb) 



<A+ 



■12a7r 2 K 2 pl 



Ttqu + avq). 



(3.25) 



The contribution from the interior of the bounce to B a is 



B, 



d x(-aR )true -Id x(-aR ) fa i se 



-2rr 2 al f d£p 3 (4:KV) 2 - I d^p 3 {AKV) 



false 



Using equations (3.12)and (3.13) we find 



B, 



_2.2 



aim 



32cwrV { / d£ 



sin 3 (i/_0 T ,2 /\, sin 3 (if + T , 2 



^ 



-Uf - rfe 



— Kl 

Hi + 



9Qan z 



kV- f kV- 9 



kV+ o / KV+ 9 

1 — r« 2+ -r^ 



(3.26) 
where V_ and V+ are the energy densities of the true and false vacua respectively. 
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Combining the contributions from the three regions, we find that to first order 



in a 



B =B - 12a7r Vpg / d(j)V2U(3U + AV ) 



+ 9Qan z 



kV- 9 / KV- o 



l _KV+-2 



, m^Pl 



(3.27) 



Now we want to determine the sign of the first order correction to the bounce 
action. If it is positive, the bubble nucleation will be less likely. First notice that 



04 



3C/V2C7 > 



(3.28) 



If the cosmological term Vq is positive, then B a \ w& \\ (3.25) has the opposite sign 
from a and 



B, 



a I wall 



.2,^2-3 



> 48|a|vr^ / pgVb / d<p V2U 



(3.29) 



Under our assumption of weak gravity, kVpq is small and we expand -B Q |i n to the 
leading order in it 



B a \ in = &<xx z K'mVo(y + -V- 



(3.30) 



This has the same sign as a if Vq is positive. Therefore the contributions to the first 
order correction from the wall region and from the bubble interior have opposite 
signs. To compare their magnitudes 



B 



a I wall 



B r 



> 



Qjd(f)V2U 
Po(V + - V-) 



(3.31) 
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notice that to leading order (see [1,3]) 

3 



po 



So we conclude 



Similarly, for 



we have 



(V+ - V.) 



"IV 



351 



(v+ - v.) 



B, 



a I wall 



B 



B r 



a 



a I in 



<0 



>2 



-B a 



, _ dB a (p ) 
dpo 



Da I wall 



B a | 

B n 



3 
> - 

2 



a 



<0 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



The bubble radius p is determined by requiring that B be stationary at that 
point. Let po be the unperturbed bubble radius and Ap the first order correction 



p = p + Ap + O(a 2 ). 



(3.36) 



We require that 



op 



= dBo 
p d P 



+ 



dB r 



Op 



To first order in a, we have 



dB 
dp 


dB 
p °P 



+ 



d 2 B n 



Pa 



dp 1 



Po 



A - _ d 2 B 
dp 2 



Ap 



po 



and 



9p 



0B a 
dp 



po 



(3.37) 



(3.38) 



(3.39) 
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Then equation (3.37) requires 

B n 



Ap 



n 



B " 



(3.40) 

Po 



As mentioned above, B a is negative when a is positive. From the work by Coleman 
and De Luccia[3] and the work by Parke [6], we know that Bq is negative. There- 
fore when a is positive, Ap is negative, which means that bubble has a smaller 
circumference radius than it does when the R 2 term is absent in the Lagrangian. 

The above conclusion can be understood intuitively. Because of the gradient 
term in the Lagrangian, the inhomogeneity in the wall region creates surface tension 
and increases the bubble energy. In order to be energetically balanced, the bubble 
interior should have a lower potential energy density to compensate the energy 
increase due to surface tension. For small a, the theory can be treated as general 
relativity with the effective potential 

V eS = V - aR 2 (3.41) 

The presence of the R 2 term lowers the potential. Both the true and false vacua 
are de Sitter spaces and the curvature scalar is proportional to the energy density. 
For the false vacuum the energy density is larger than for the true vacuum, so the 
relative energy difference is reduced by the presence of the R 2 term. However in 
the thin wall region the curvature scalar becomes fairly big when the field changes 
rapidly. This means that in the wall region the potential energy decreases at a 
much larger rate than the potential energy inside the bubble. So both the surface 
tension and the interior potential energy difference decrease, but the surface tension 
decreases at a larger rate. So to be balanced energetically, the volume inside the 
bubble should decrease. This is why the bubble radius decreases. Generally, it 
is easier for a small bubble to nucleate, so the bounce action decreases and the 
bubble nucleation rate increases. 
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4. Bounce Action: Conformal Transformation Method 

In this section, we are going to use the conformal transformation of Section II 
to study the bubble nucleation rate in fourth-order gravity. If the Lorentzian action 
is Eq. (2.17), then the Euclidean action is 



S E = / <ftx\[g 



i^- rf+ >« !+,, («i' 



(4.1) 



The conformal transformation described in Section II leads to the action 



S 



E 



J4 ~ 

a xy g 



m: 



167T 



-R 



Ylo?Kd ix ad^a d^<p d»<p V{4>) + aa 2 



[1 + Aetna) 2 2(1 + Aetna) (1 + Aetna) 2 



(4.2) 

Assuming 0(4) symmetry, we use the metric (3.4). The Euclidean action then 
becomes 



S E = 2vr 2 I ddp 3 



2(1 + Aetna) 

/2 



+ 



12ct 2 na' 2 V + eta 2 



3 / 2 // i '^ 
+ -\P P + PP ~ P) 



K 



1+ Aetna) 2 (1+ Aetna) 2 



(4.3) 



The Einstein equation is 



/2 , . 1 2 
P = 1+ g^P 



/,/2 



12a 2 Ka /2 



aa 



y 



2(l + 4aKcr) (1 + Aetna) 2 (1 + 4a«;cr) 2 



(4.4) 



The other equations of motion are 



i" + £^0' = ^ : 

p 1 + 4a«a 



(4.5) 



and 



,/2 jl/2 



1 a - 4kV 



„ 3p' , 8o;Kcr'' 

a" + -i-cr' = — + 

p 1 + Aetna 6a 6a/t 1 + 4aK<7 

When ex. — 0, these equations of motion agree with Eqs. (3.5) and (3.6). 



(4.6) 
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The second derivative of p can be eliminated by integration by parts. (The 
surface term from the parts integration is harmless because we are only interested 
in the action difference between two solutions that agree at the boundary.) We 
thus obtain 



S E = 2tt 2 I ddp 3 



1/2 



12a 2 Ka f " 



--(PP' 2 + P)\- 

K 



V + aa 2 



2(l + 4a«tr) (l + Aana) 2 (1 + Aetna) 2 



(4.7) 



We now use Eq. (4.4) to eliminate p'. We find 



S E = 4vr z \ di 



3 3 V + aa 2 

k (1 + Aana) 2 



(4- 



Now we evaluate the above integral in three regions. Outside the bubble the 
bounce solution agrees with the false vacuum and the contribution to the bounce 
action is zero: 



B, 



out 







(4.9) 



The interior of the bubble is de Sitter space and from the equations of motion we 
know 



a = 4kV 
sin(#0 



(4.10) 



P 



H 



with 



iT 



kV 



3(l + 16a/t 2 \/)' 



(4.11) 



The contribution to the bounce action is the difference of the integral (4.8) between 
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the true vacuum and false vacuum. It is easy to find that 
12tt 2 



An = -^r(l + l&wcV) { [I - (Hpff 2 - l) - (/aise) (4.12) 

K, Z V I L J J true 



Expanding this to first order in a, we find 

Bin =Bf Q (p) 



+ 96an 2 



kV- „ / kV- o 
1 - ~^P 2 2 + — - p 2 



kV+ f kV+ _ 2 

1 — 7T P 2 ^ — ^~P 



(4.13) 



which is exactly the same result as in the previous section. 

In the wall region, we expand the integral (4.8) to first order in a. This gives 

5wall = 4vr 2 p 3 j di{V b - V f ) + 47r 2 ap 3 J d£ [(a 2 - 8KaV) bounce - (false)} (4.14) 

At first sight one might think that the first integral is the unperturbed contribution 
with a = and the second integral is the first order correction. However, the first 
integral also contains a first order correction because V is a function of <fi and when 
a^O the dependence of on £ changes. Let us investigate the equations of motion 
as a power series in a. From Eqs. (4.6) and (4.5), the equations of motion in the 
thin-wall approximation are 



a = k{4V + <f)' ) + 0{a) 



dV 



+ Aana 



dV 



Aetna—— + 0(a z 



(4.15) 



Multiplying the second equation by <p', we get 



V = AOLK 

To first order in a, this implies 

!,/2 



a'6' 2 - aV 



V - V + 0(a) 



0(a z 



(4.16) 



(4.17) 
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Together with the first equation in (4.15), this leads to 

a'<j)' 2 - aV' = k(3U 2 - AVqU)' + 0(a) (4.18) 

Putting this back in (4.16) and integrating, we get 

<j)' 2 = 2U(1 + 12aK 2 U - lQaK 2 V ) + 0(a 2 ) (4.19) 

The first term of (4.14) is 

4vr 2 p 3 J dt (H - V f ) = 4vr 2 p 3 J^t(V b -V f ) (4.20) 

Using Eq. (4.19) , we find the first order correction to this is 

4tt 2 p 3 / d<j) aK 2 V2U(-3U + AVq) (4.21) 



The contribution to the first order correction from the second integral is easily 
found to be 

47rVW I d(j> V2U(-6U - 16Vb) 
Therefore the first order correction due to the wall region is given by 

5 Q |wall = -127c 2 aK 2 p 3 J dcf) V2U(3U + AV ) (4.22) 

Again, this is exactly the same result as in last Section. Although p is determined 
in the new metric now, the difference between the old metric and the new one is 
first order in a. So it does not invalidate our conclusion. 

Therefore, we get exactly the same results by using two completely different 
methods. This also proves the validity of the conformal transformation we dis- 
cussed in Section II. 

I thank Erick Weinberg for suggesting this subject and for constant discussion 
and encouragement. 
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